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Abstract  By  using  different  isospin  dependent  proximity-type  potentials,  we  performed  a 
detailed  study  of  neutron/proton-rich  colliding  nuclei  with  N/Z  ratio  between  0.5  and  2.0. 
Isotopes  of  three  different  series  namely,  Ne-Ne,  Ca-Ca,  and  Zr-Zr  are  taken  into  account.  A 
monotonous  increase  (decrease)  in  the  fusion  barrier  positions  (heights)  using  a  unified  second 
order  nonlinear  parametrization  in  the  normalized  fusion  barrier  positions  and  heights  with 

—  1  j  is  presented.  These  predictions  are  in  good  agreement  with  the  available  theoretical 

as  well  as  experimental  results.  The  fusion  probabilities,  however  shows  linear  dependence. 
Further,  the  neutron/proton  content  plays  dominant  role  at  near  barrier  energies  only.  Our 
results  are  also  independent  of  the  model  used  as  well  as  reaction  partner  and  isospin  content. 
More  experiments  are  needed  to  verify  our  predictions. 
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1.  INTRODUCTION 

With  the  availability  of  radioactive-ion  nuclear  beams  the  fusion  of  colliding  nuclei  with 
the  excess  of  neutron/proton  ratio  or  near  the  drip  line  has  attracted  central  position  in  the 
current  days  research  [1-14].  The  neutron-rich  radioactive-ion  beams  have  been  applied 
to  synthesize  new,  neutron-rich  heavy  nuclei  [1].  This  is  because  in  the  synthesis  of  heavy 
nuclei  with  n-rich  projectile,  one  expect  a  higher  survival  probability  of  the  completely 
fused  system  due  to  its  lower  hssility  and  lower  excitation  energies.  One  of  the  important 
use  of  radioactive-ion  beam  is  to  synthesize  and  study  super  heavy  elements  [2,  6]. 

Several  experiential  studies  have  also  been  undertaken  in  the  literature  to  study  the 
effect  of  varying  the  N/Z-ratio  of  the  projectile  and  target  nuclei  upon  the  fusion 
cross  sections  [7-10,  15-21],  Recently,  the  enhancement  of  the  fusion  cross  section 
has  been  observed  for  neutron-rich  projectile  [18-21].  This  is  due  to  the  fact  that 
neutron  density  distribution  extended  far  beyond  the  range  of  normal  nuclei  leading 
to  a  larger  barrier  radius  and  hence  to  a  lower  height  of  the  fusion  barrier.  Similarly, 
Bian  et  al.  [7],  employed  the  concept  of  static  and  dynamical  Coulomb  barrier  to 
analyze  the  enhancement  of  the  fusion  cross  section  for  neutron-rich  systems  using 
IQMD  model.  It  has  been  suggested  that  the  fusion  yield  would  be  further  enhanced 
when  the  reaction  is  induced  by  unstable  neutron-rich  nuclei  [22],  The  neutron 
flow  effect  were  also  not  proposed  to  understand  the  observed  enhancement  of  the 
fusion  cross  sections  [23,  24].  In  addition,  the  fusion  transfer  at  the  neck  region 
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has  been  suggested  [25,  26].  A  microscopic  description  of  the  formation  of  neck 
in  the  fusion  reaction  remain  a  challenge  to  microscopic  theories.  The  influence 
on  sub-barrier  fusion  of  processes  such  as,  transfer  [27]  and  breakup  reaction 
[14,  28]  is  not  yet  clear;  moreover,  the  effect  of  unusual  structure,  such  as  halos  and 
skins  [29],  is  bing  studied  [8,  30].  Recently,  Sun  etal.  [31],  suggested  that  N/Z  may  be 
used  as  an  experimental  observable  to  extract  neutron  skin  for  neutron  rich  nuclei.  All 
above  experimental  as  well  as  theoretical  information  indicates  that  the  dynamics  of 
neutron/  proton-rich  nuclei  is  not  fully  understood  and  needs  further  attention. 

The  properties  of  various  neutron-rich  nuclei  with  different  N/Z  -ratio  are  studied  in 
the  literature  e.g.:  9~102He( N/Z  =3.50-4.00;  where  N  and  Z  are  the  neutron  and  proton 
content  of  the  nucleus),  6-8'9'u3Zi  (N/Z  =1.0,  1.67,  2.0,  2.67), 22 f  (N/Z  =2.67),  2("2\0 
(N/Z  =  2.25-2.50),  2\F  (N/Z  =  2.444), 28- 3410/V<?  (N/Z  =1.8,  2.2,  2.40),  30-32'37u/Va  (N/Z 
=  1.727-1.909,  2.364),  40nMg  (N/Z  =  2.333),  49-51lgAr  (N/Z  =  1.722-1.833),  6020Ca  (N/Z  = 
2.0 ),57'6025Mn  (N/Z  =  1.28-1.4),  68-782gM  (N/Z  =  1.429-1.786),  84'8630Zh  (N/Z  =  1.8,  1.87), 
9a  9232Ge  (N/Z  =  1.813,  1.875),  m5(jSn  (N/Z  =  1.64),  12347Ag  (N/Z  =  1.617),  m'm4gCd 
(N/Z  =  1.563-1.667)  [32]  and  proton-rich  are 6 fie  (N/Z  =  0.50), 10  (N/Z  =  0.429),  1280 
(N/Z  =  0.50),  179F  (N/Z  =  0.89),  2214Si  (N/Z  =  0.571),  31lgAr  (N/Z  =  0.722),  \Ca  (N/Z 
=  0.70),  38'392,77  (N/Z  =  0.727,  0.773),  45'492/e  (N/Z  =  0.731,  0.885),48’49'532gM  (N/Z  = 
0.714-0.75,  0~893),  5430Zn  (N/Z  =  0.80),  217g2U  (N/Z  =  1.359)  etc.  [33,  34], 

A  suitable  set  of  models  are  therefore  needed  to  study  the  dynamics  of  neutron/  proton- 
rich  colliding  nuclei.  A  large  number  of  theoretical  models  are  available  in  the  literature 
based  upon  the  different  assumptions  [2-6].  Among  them,  proximity  potential  due  to 
Blocki  et  al.  [35],  is  well  known  for  its  simplicity  and  wider  applications  in  different 
fields.  Several  modifications  or  refinements  over  the  original  proximity  potential  are 
also  available  in  the  recent  time  by  including  either  up-to-date  knowledge  of  the  surface 
energy  coefficient  or  nuclear  radii  [4,  36].  Various  authors,  modified  or  parametrized 
their  approaches  within  the  proximity  concept  [4].  All  these  modifications  include  new 
emerging  degree  of  freedom  i.e.  isospin,  either  in  radius  formula,  universal  function 
and/or  in  surface  energy  coefficient  [4].  Further,  The  outcome  will  definitely  different 
if  one  use  such  type  of  models  in  the  isospin  plane. 

Recently,  one  of  us  and  collaborator,  have  carried  out  a  detailed  study  involving  large 
number  of  symmetric  as  well  as  asymmetric  colliding  nuclei  using  16  proximity-type 
potentials  [4],  Unfortunately,  the  maximum  N/Z  content  of  all  the  experimentally 
studied  heavy-ion  reactions  is  1.60  (i.e.  6He  +  238 //)  [8].  On  the  other  hand,  the  first 
measurement  with  the  proton  drip  line  nucleus  is  of  11 F  +  2S2Pb  (with  N/Z  =1.473)  [33]. 
Therefore,  a  systematic  study  of  nuclei  having  larger  N/Z  ratio  using  new  proximity- 
type  potentials  is  in  demand.  Further,  it  gives  us  a  unique  possibility  to  test  the  validity 
or  accuracy  of  these  models  for  the  nuclei  far  away  from  the  line  of  stability.  Therefore, 
a  systematic  dependence  of  fusion  barrier  (heights  and  positions)  and  cross  sections 
using  various  isospin  dependent  models  on  neutron  excess  is  needed.  Similar  study  was 
also  presented  by  Puri  et  al.  [3],  where  only  Ca  and  Ni  series  were  used.  In  the  present 
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study,  we  extend  the  work  to  include  new  series  like  Ne-Ne  (with  N/Z  ratio  =  0.6-2.0) 
and  Zr-Zr  (with  N/Z  =  0.75-2.0)  along  with  Ca-Ca  (with  N/Z  =  0. 5-2.0)  series.  The 
overall  domain  of  N/Z  ratio  is  from  0.5  to  2.0  for  all  series.  The  asymmetry  parameter 
A  (=(MZ)-1)  of  the  colliding  nuclei  varies  between  -0.5  and  1.0.  Note  that  non  zero 
value  of  A  will  involve  complex  interplay  of  the  isospin  degree  of  freedom  which  has 
strong  role  at  intermediate  energies  as  well  [37].  Section  2,  deals  with  fine  points  of  the 
models  in  brief.  Section  3  contains  the  results  and  summary  is  presented  in  Section  4. 
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2.  FORMALISM 


The  total  ion-ion  interaction  potential  V  (r)  comprises  of  nuclear  and  Coulomb  part: 


VT(r)  =  VN(r)+Vc(r).  (1) 

Here  Vc(r)  =  (ZjZ2e2  /  r)  is  a  good  approximation,  because  fusion  happens  at  a  distance 
greater  than  touching  configuration  of  the  colliding  pair. 


The  nuclear  part  of  the  ion-ion  interaction  potential  V  (r)  is  calculated  within  the 
proximity  concept.  All  proximity  potentials  are  based  upon  the  proximity  force 
theorem,  according  to  which  [35],  “ the  force  between  two  gently  curved  surfaces 
in  close  proximity  is  proportional  to  the  interaction  potential  per  unit  area  between 
the  two  flat  surfaces”.  In  original  proximity  potential  [35],  the  nuclear  part  of  the 


interaction  potential  V  ( r)  can  be  written  as 


VN(r)  =  Ait^bCQ 


r  —  C,  —  C, 


MeV. 


In  this,  C  is  reduces  radius  with  equivalent  sharp  radius  R  as 


(2) 


R.  =  1.28 A‘n  -  0.76  +  0.8 A/1'3  fm  (i  =1,2).  (3) 

r—C  —C 

Where  $(£  = — 'h  2 )  is  the  universal  function  that  depends  on  the  separation 
between  the  surfaces  of  two  colliding  nuclei  only.  Both  these  factor  do  not  depend 
on  the  isospin  content.  However,  the  last  parameter  7,  the  surface  energy  coefficient, 
depends  upon  the  neutron/proton  excess  as 


7  =  7 


0 


N-Z) 
A 


(4) 


where  N,  Z  being  the  total  number  of  neutrons  and  protons.  In  the  original  version,  70 
=  0.9517  MeV/fm2  and  ks  =  1.7826  [35].  Noted  that  for  symmetric  nuclear  matter,  N  = 
Z,  7  =  70  =  0.95 17  MeV/fm2  indicating  maximum  strength  of  the  potential.  If  we  move 
to  neutron  (proton)-rich  colliding  nuclei  with  N  >  Z  (N  <  Z)  then  7  starts  decreasing 
resulting  in  comparatively  lesser  attractive  potential. 


Later  on,  these  coefficients  were  further  improved  by  Moller  and  Nix  with  values  70 
=  1.2496  MeV/fm2  and  ks  =  2.3  [38].  This  is  labeled  as  Prox  88  [4,  39].  In  the  latest 
version  of  proximity  potential  [36],  7  has  form  based  on  the  precise  neutron  skin  as 
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7  =  ' 


47T  r„ 


18.63(MeV)-Q 


2r; 


(5) 


The  corresponding  proximity  potential  is  labeled  as  Prox  00  [4],  One  of  us  and 
collaborator  [4],  modified  above  potential  to  include  latest  radius  formula  [40]  and  is 
denoted  as  Prox  00DP.  Note  that  both  Prox  00  and  Prox  00DP  has  isospin  dependent 
radius  with  slightly  diffierent  constants  whereas  the  factors  surface  energy  coefficient 
7  and  the  universal  function  <h(.y)  are  same.  In  both  newer  versions  of  Bass  (labeled  as 
Bass  77  and  Bass  80  in  Ref.  [4])  radius  is  slightly  changed  to 

R.  =  1 .  16A1/3  - 1 .39A71'3  fm  (i  =  1, 2),  (6) 

and  then  to  sharp  radius  as  is  used  in  Prox  77  in  the  later  version  (i.e.  Bass  80).  Both 
newer  versions  of  Winther  (marked  as  BW  91  and  AW  95  in  Ref.  [4])  has  again  similar 
expression  for  7  as  is  given  in  Prox  77  with  slight  diffirence  that  here  isospin  content 
is  calculated  separately  for  the  target/  projectile.  Whereas,  first  version  due  to  Winther 
(labeled  as  CW  76  in  Ref.  [4])  does  not  have  any  7  dependence.  Even  radii  are  function 
of  mass  only.  Both  versions  of  Ngo  (labeled  as  Ngo  75  and  Ngo  80  in  Ref.  [4])  do 
not  consider  7,  but  latest  version  of  Ngo  (Ngo  80)  has  isospin  dependence  in  radius 
parameter.  On  the  other  hand,  a  complex  isospin  dependence  in  the  universal  function 
$(s)  and  radius  is  given  in  the  version  of  Denisov  [41  ].  Also  by  using  the  latest  form 
of  radius  given  in  Ref.  [40]  in  Denisov  potential  resulting  in  closer  agreement  with 
the  experimental  data  for  fusion  barrier  heights  and  cross-sections.  This  potential  is 
labeled  as  Denisov  DP  [4] .  All  the  above  mentioned  proximity-type  potentials  are  able 
to  reproduce  the  experimental  fusion  barriers  within  ±10%  on  the  average  [4], 


In  total,  eight  proximity-type  potentials  are  used  in  the  present  study.  Among  them, 
three  are  basic  proximity  potentials  (Prox  77,  Prox  88,  and  Prox  00),  Bass  potential 
(Bass  80),  Winther  (AW  95)  and  Ngo  (Ngo  80)  each,  and  two  newly  modified 
potentials  (Prox  00DP  and  Denisov  DP)  are  used.  The  model  due  to  Bass  et  al., 
(Bass  80)  used  in  the  present  analysis  is  independent  of  isospin  dependence.  For  the 
detail  of  the  models  reader  is  refer  to  Ref.  [4],  In  most  of  the  potentials  and  versions, 
modifications  are  made  either  through  the  surface  energy  coefficients  or  in  nuclear 
radii.  These  two  rather  being  technical  parameters  can  have  sizeable  effects  on  the 
outcome  of  a  reaction  [5]. 


From  these  brief  outlines,  it  is  clear  that  much  stress  is  made  on  the  surface  energy 
coefficients  7  and  nuclear  radii  to  incorporate  the  isospin  factor  of  a  potential.  Definitely, 
based  on  different  assumptions  and  isospin  dependence,  different  versions  will  respond  to 
the  collision  of  neutron-rich  or  -deficient  nuclei  differently  compared  to  N  =  Z  nuclei. 

3.  RESULTS  AND  DISCUSSIONS 

The  present  study  deals  with  large  variety  of  above  mentioned  potentials.  Using  these 
potentials,  we  firstly  calculate  the  total  ion-ion  interaction  potential  using  Eq.  (1). 


162 


Once  total  ion-ion  interaction  potential  is  calculated,  one  can  extract  the  barrier  height 
VB  and  barrier  position  Rg  using  conditions: 


dVT(r)  d2V  (r ) 

„  =  0;  and 

•r—Ko 


R  <0. 

-Kd 


(7) 


Here 
A 


dr  8  dr 

we  consider  the  collisions  of  three  different  series  namely; 
'Ca  +  ^ Ca  (with  Nz  =  0.5  to  2.0) ,  and  A'Zr  +  AlZr  (with  N?  =  0.75  to  2.0)  to  cover 
wider  mass  range.  We  starts  with  the  collision  of  N  =  Z  nuclei  and  then  add  (or  remove) 

neutrons  gradually  from  either  of  the  colliding  pairs  till  we  reach  N  =  2Z  (N  =  0.5Z) 
nuclei.  In  total,  150  such  collisions  involving  different  isotopes  of  different  series  are 
taken  into  account. 
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As  a  first  step,  we  check  the  effect  of  addition  or  removal  of  neutrons  on  the 
nuclear  part  of  the  interaction  potentially  (r)  in  different  models.  In  Fig.  1, 
we  display  V  (r)  as  a  function  of  internuclear  distance  r  for  the  reactions  of  70Zr  + 
70Zr,  ™Zr  +  76Zr,  16Zr  +  mZr,  mZr  +  mZr,  mZr  +  S6Zr,  S6Zr  +  90Zr,  and  90Zr  +  90Zr  using 
eight  proximity-type  potentials.  Based  upon  the  different  assumptions  used  in  different 
models,  shape  as  well  as  strength  of  the  potential  differ  accordingly.  At  the  same  time, 
the  general  shape  at  the  surface  region  is  same.  In  particular,  Bass  80,  Prox  00,  and  Prox 
00DP  have  no  repulsive  core  at  shorter  distances,  whereas  AW  95  follow  Woods-Saxon- 
type  form.  On  the  other  hand,  Ngo  80,  Denisov  DP,  Prox  77,  and  Prox  88  have  repulsive 
core  at  shorter  distances.  In  addition,  four  effects  of  addition  (removal)  of  neutrons  to 
N  =  Z  nuclei  are  clearly  visible:  (i)  barrier  height  is  decreased  (ii)  barrier  position 
is  increased  (iii)  depth  of  the  potential  is  increased  in  all  potentials  except  Denisov 
DP,  and  (iv)  diffuseness  of  the  potential  is  also  changes.  These  effects  will  definitely 
influence  the  fusion  probability  at  below  barrier  energies. 

Therefore,  before  discussing  the  enhancement  of  the  fusion  cross  section  for  the 
neutron-rich  fusion  reaction,  we  investigate  the  systematic  dependence  of  fusion 
barriers  on  the  neutron  asymmetry  parameter  Ay  =  £  —  lj. 

Using  all  above  stated  models,  fusion  barrier  heights  and  positions  are  calculated 
for  more  than  150  colliding  pairs.  The  variation  in  the  fusion  barrier  positions  and 
heights  with  neutron/proton  content  is  analyzed  in  Figs.  2  and  3,  respectively.  Here, 
we  display  the  variation  of  A RB  (%)  and  AVg  (%)  defined  as; 

^ 

AR„(%)  =  — - -  x  100, 

8  K 

v  -y° 

AV_  (%)  =  — - -  x  100, 

B  y° 

T  » 


(8) 

(9) 
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Figure  1:  The  nuclear  part  of  the  interaction  potential  V  (MeV)  as  a 
function  of  internuclear  distance  r  for  the  reactions  of  70Zr  +  ™Zr,  70Zr  + 
76Zr,  76Zr  +  mZr,  mZr  +  mZr,  mZr  +  86Zr,  86Zr  +  90Zr,  and  90Zr  +  90Zr  using 
eight  sets  of  proximity-type  potentials. 


against  asymmetry  parameter  As  (=  A  ~  l)  using  eight  set  of  potentials  discussed  above. 

Here,  R°  and  VB  are,  respectively,  the  fusion  barrier  position  and  height  corresponding  to 
(N  =  Z)  colliding  nuclei  and  RB  and  VB  refer  for  neutron/proton-rich  colliding  nuclei.  The 
main  advantage  of  these  normalized  variation  is  that  it  gives  mass  independent  picture.  For 
the  present  picture,  we  starts  with  the  collision  of  N  =  Z  nuclei,  then  gradually  add/remove 
neutrons  from  either  of  the  colliding  pair.  For  example,  we  started  with  the  collision  of 
mCa  +  40Ca,  then  add  neutrons  gradually,  by  keeping  charges  Zl  and  Z,  always  fixed.  In  this 
series,  at  the  end  of  the  chain  we  have  the  collision  of  mCa  +  60  Ca.  Similarly,  if  we  remove 
the  neutrons  then  we  have  at  the  end  the  collision  of  30Ca  +]llCa.  It  is  clear  from  the  figures, 
that  all  the  models  follow  a  unified  non  linear  second  order  parametrization  given  as: 


A  R. 


AV„ 


N 

fN 

(%)  =  a 

- 1 

+  b 

- 1 

Z 

KZ 

(N 

(N 

(%)  =  c 

- 1 

4-  d 

- 1 

Z 

Z 

(10) 

(11) 
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reverse  happen  for  neutron  -deficient  nuclei.  All  models  do  not  show  much  scattering  from 
the  middle  curve  as  one  move  from  N  =  Z  to  very  neutron-rich  and  -deficient  colliding 
nuclei.  Bass  80  and  Ngo  80  show  slight  scattering  for  neutron-rich  and  -deficient  nuclei.  It 
may  be  due  to  the  reason  that  the  isospin  dependence  included  in  Ngo  80  model  have  not 
much  effect  on  [ :V  j  ratio,  whereas,  Bass  80  is  independent  of  such  kind  of  dependence.  We 

\  Z) 

further  note  that  the  slope  of  the  central  lines  also  varies  from  model  to  model,  whereas, 
the  overall  pattern  is  nicely  explained  by  the  above  parametrizations.  It  may  be  due  to  the 
different  assumptions  used  in  different  models. 


Here,  a,  b,  c,  and  d  are  the  constants  varies  from  model  to  model  and  its  values  are  displayed 
in  Figs .  2  and  3 .  The  above  results  are  in  good  agreement  with  the  work  due  to  Puri  el  al. ,  for 
Ca  and  Ni  series  [3] .  The  available  experimental  as  well  as  theoretical  data  is  also  displayed. 
All  the  available  theoretical  as  well  as  experimental  data  follow  our  parametrization  pattern 
very  well  except  few  points  due  to  experimental  uncertainty  in  different  experimental 
setups.  The  theoretical  data  is  taken  from  Refs.  [42],  whereas,  the  experimental  data  is 
taken  from  Refs.  [43].  The  above  pattern  indicates  that,  with  the  addition  of  neutron  to 
N  =  Z  nuclei,  fusion  barrier  position  is  increased  and  barrier  height  decreased,  whereas. 


Along  with  A RB  (%)  and  AVg  (%),  we  have  also  studied  AVv  (%)  and  AFC  (%)  using 
the  same  set  of  models  and  series  (not  shown  here).  We  see  that  these  variations  show 
larger  scattering  due  to  the  different  assumptions  resulting  from  different  forms  of 
radii,  universal  function,  surface  energy  coefficients  etc.  The  diffuseness  parameter 
is  also  different  in  some  potentials. 


We  further  analyze  the  effect  of  addition  or  removal  of  neutrons  on  fusion 
probabilities.  In  Fig.  4,  we  display  Act(us  (%)  with  As  L--ij  at  energy  E  =  1.02  VB°. 
Where  A<rfus  (%)  is  calculated  similar  to  Eqs.  (8)  and  (10)  as 


a,  (E°  )-a°  (E°  ) 
A(7fJ%)  =  J  x  10Q> 

<Uus(£cJ 


(12) 


with  £c°m  and  <r°s  being  the  center-of-mass  energy  and  fusion  cross  section  for 
symmetric  colliding  pair.  From  this  figure,  it  is  clearly  understood  that  all  models 
follow  linear  dependence  in  fusion  probabilities  with  As  of  the  form. 


Aabis%  =  a 


(13) 


Here  a  is  a  constant  and  varies  from  model  to  model.  Its  values  are  displayed  in  the 
figure.  The  overall  variation  in  its  values  is  from  721.58  ±  40  to  1036.11  ±  150.  It 
indicates  that  some  potentials  strongly  depend  upon  the  N/Z  content  of  the  reaction. 
The  shaded  part  represent  the  slight  deviations  from  the  linear  dependence.  As  we 
noticed  earlier,  the  slope  of  the  linear  curve  is  changes  as  one  moves  from  one  model 
to  other.  But  the  general  behavior  is  nicely  explained  by  linear  parametrization  [Eq. 
(13)].  It  clearly  indicates  that  our  results  are  system  as  well  as  model  independent. 
Also,  the  isospin  degree  of  freedom  introduced  in  different  potentials  does  not  affect 
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Figure  2:  The  normalized  barrier  positions  A RB  (%)  [Define  in  Eq.  (8)]  as  a 
function  of  As  {  =  A  —  1  j  We  display  the  results  of  our  calculations  for  the 
collisions  of  Ai  Ne  +  Al  Ne, A>  Ca  +  Ca,  and  A<  Zr  +  A2  Zr  series  using  eight 
models  along  with  other  available  theoretical  and  experimental  values.  The 
theoretical  as  well  as  experimental  data  reported  here  is  taken  from  Refs. 
[42]  and  [43],  respectively.  The  shaded  areas  denotes  the  deviation  from  the 
cental  solid  lines. 


our  outcome.  We  carried  out  similar  study  at  different  energies  i.e,  E  =  1.075  VB,  1.125 
VB°,  1.25  VB°,  1.375  VB°,  and  1.50  VB°.  We  see  that  the  effect  at  near  barrier  energy 
is  more  pronounced,  whereas,  at  above  barrier  energies  the  effects  are  insignificant. 
We  further  test  our  parametrized  results  [Eqs.  (10)  and  (11)]  in  Fig.  5  by  using  the 
formula  due  to  Puri  et  al.  [3]  for  neutron/proton-rich  colliding  pair  as 


<,s(mb)  =  107rR°B(l  +  ARB) 


1- 


VB°(1  +  AVB) 


(14) 
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Figure  3:  Same  as  Figure  2,  but  for  AKfi  (%)  [Define  in  Eq.  (9)]. 


for  Ca-Ca  series  only.  The  corresponding  results  obtained  using  sharp  cutoff  formula 


<7fuS(mb)  =  1°7rRB 


(15) 


are  also  displayed.  From  the  figure,  it  is  clear  that  our  parametrized  results  are  very  close 
to  the  actual  one  obtained  using  Eq.  (15)  and  is  in  agreement  with  earlier  predictions  [3J. 
This  is  the  first  time  that  such  a  wider  mass  range  is  taken  into  account  and  a  complete 
isotopic  dependence  is  studied  using  different  proximity-type  potentials.  This  study  will 
certainly  helpful  to  study  various  neutron/proton-rich  nuclei  lying  between  these  limits. 


4.  SUMMARY 

We  analyzed  the  fusion  of  three  different  series  namely,  Ne,  Ca,  and  Zr  by  covering 
the  wider  mass  spectrum  with  |^)  ratio  between  0.5  and  2.0.  We  analyzed  the 
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Figure  4:  The  normalized  fusion  probabilities  A(jfus  (%)  as  a  function  of 
As[=^  —  l)  for  the  collisions  of -hNe  +  ^Ne,  AlCa  +  A*Ca,  and  A'Zr  +  AiZr 
series  using  eight  sets  of  proximity-type  model  at  incident  energy  Ec  =  1 .02 

V* 

systematic  dependence  of  fusion  barriers  on  neutron  excess  and  presented  a  unified 
second  order  non  linear  quadratic  parametrization  in  fusion  barrier  heights  and 
positions  with  (-  —  l)  using  eight  isospin  dependent  proximity-type  models  for  three 
different  series.  Our  results  are  in  good  agreement  with  the  available  theoretical  as 
well  as  experimental  results.  A  linear  dependence  in  the  fusion  probabilities  is  also 
presented.  Further,  the  enhancement  in  fusion  cross  sections  for  neutron-rich  nuclei 
due  to  lowering  of  fusion  barrier  heights  is  clearly  seen,  whereas,  reverse  happen  for 
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Figure  5:  The  parametrized  fusion  cross  sections  <rfus(nib)  (symbol)  using 
[Eq.  (14)]  and  the  exact  calculations  (line)  using  [Eq.  (15)]  verses  the 
center-of-mass  energy  E  for  the  reactions  of  uCa  +  MCa,  uCa  +  i0Ca, 
mCa  +  40Ca,  40Ca  +  48Ca,  and  48Ca  +  48Ca  using  eight  sets  of  proximity- 
type  potentials. 

proton-rich  nuclei.  Along  with  this,  our  parametrization  pattern  is  independent  of 
the  colliding  nuclei  as  well  as  model  and  isospin  content.  At  near  barrier  energies, 
N/Z  content  plays  dominant  role,  whereas,  the  effect  is  insignificant  at  above  barrier 
energies.  More  experiments  are  needed  to  verify  our  predictions. 
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